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Abstract 



We establish the existence and uniqueness of both local martingale and local pathwise solutions of an 
abstract nonlinear stochastic evolution system. The primary application of this abstract framework is 
2 ' to infer the local existence of strong, pathwise solutions to the 3-D primitive equations of the oceans and 

atmosphere forced by a nonlinear multiplicative white noise. Instead of developing our results specifically 
for the 3D primitive equations we choose to develop them in a slightly abstract framework which covers 
many related forms of these equations (atmosphere, oceans, coupled atmosphere-ocean, on the sphere, 
on the /3-plane approximation etc and the incompressible Navier-Stokes equations). In applications, all 
of the details are given for the /3-plane approximation of the oceans equations. 
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CN ■ 1 Introduction 

(•~^ ■ In this work we develop a local existence theory for a class of abstract stochastic evolution systems of the 

f^ , form: 

dU + {AU + B{U,U)+F{U))dt^a{U)dW, U{0)^Uo, (1.1) 



which includes the primitive equations of the ocean as explained below. We choose to develop our setting 
in a somewhat abstract framework so that our results may cover some closely related equations such as the 
primitive equations of the atmosphere or the coupled atmospheric/oceanic system, equations with chemistry, 
equations on the sphere or /3-plane approximations, etc. We do not expand further on these latter applications 
in this article in order to avoid excessive developments. The stochastic primitive equations of the ocean have 
been previously studied in [26, 22, 24, 15, 27] but none of these works address the full 3-d system in the 
context of a nonlinear multiplicative noise. 

The deterministic primitive equations are widely seen as a fundamental model for large scale oceanic and 
atmospheric systems. For the oceans they are derived from the fully compressible Navier-Stokes equations 
combined with the Boussinesq and hydrostatic approximations. See e.g. [39] for further physical background. 
Given the growing importance of probabilistic methods in sub-grid scale parameterization, geophysicists 
currently devote significant attention to stochastic forms of the equations of geophysical fiuid dynamics; see 
e.g. [14, 40, 41, 44, 32, 34, 4, 50]. From the mathematical point of view the stochastic primitive equations 
have been considered for a two dimensional version of the equations in [26, 15, 25]. In three space dimensions 
[27] has addressed the case of an additive noise. In this situation the outcome w (in the underlying probability 
space fi) may be treated as a parameter in the problem. Furthermore the methods in [27] do not allow for 
physically realistic boundary conditions which we are able to treat here with our methods. 
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Let us recall at this point that the theory of the related stochastic Navier-Stokes equations have undergone 
substantial developments; see e.g. [3, 48, 9, 8, 17, 12, 2, 5, 6, 19, 10, 36, 37, 45, 18, 23]. However to 
emphasize the notable differences between the stochastic Navier-Stokes Equations and the primitive equations 
we recall that the deterministic Primitive equations arc known to be well posed in space dimension three 
(see [7, 30, 31]). This result is not known for the Navier-Stokes Equations and is the object of the famous 
Clay problem. On the other hand the primitive equations are technically more involved than the Navier- 
Stokes equations. For further background concerning the mathematical theory for the deterministic primitive 
equations see e.g. the review article [42] and the references therein. The stochastic primitive equations that 
wc consider in this work are described in detail in Section 6 below. 

In the theory of stochastic evolution equations two notions of solutions are typically considered namely 
pathwise (or strong) solutions and martingale (or weak) solutions. In the former notion the driving noise 
is fixed in advance while in the later case these underlying stochastic elements enter as an unknown in the 
problem. In this work we will consider both notions for (1.1) and illuminate the relationship between these 
two types of solutions. The classical Yamada-Watanabe theorem from finite dimensional stochastic analysis 
says that pathwise solutions exist whenever martingale solutions may be found and pathwise uniqueness 
of solutions holds (see e.g. [43]). Similar results have later been established along more elementary lines, 
in [28], using a simple characterization of convergence in probability (see Proposition 2.2 below). This 
characterization may also be employed in the infinite dimensional context and is used below to pass from 
the case of martingale to pathwise solutions. In any case to the best of our knowledge, no one has previously 
established such a 'Yamada-Watanabe' type result for the Primitive Equations (or for that matter for the 
3-D Navier-Stokes equations). 

The exposition is organized as follows. In Section 2 we make precise the set-up of the abstract problem 
(1.1) and briefly recall some relevant mathematical preliminaries from probability theory and functional 
analysis. A Galerkin scheme for (1.1) is considered in Section 3. By making use of an appropriate cut- 
off function in the formulation of the equations we are able to establish uniform a-priori estimates for the 
corresponding sequence of approximate solutions. In Section 4 we outline the compactness arguments that 
lead to the local existence of martingale solutions. We turn then to pathwise solutions in Section 5. Here 
the first step is to establish conditions for pathwise uniqueness. We then revisit the compactness methods 
described in the previous section now making use of the additional criteria for convergence in probability. 
In Section 6 we apply the abstract results to the stochastic primitive equations. In the final Sections, 7, 8, 
we provide the technical details of the passage to the limit from compactness of Galerkin approximations 
established earlier in Sections 4, 5. 

2 Mathematical Framework 

In this section we set up the abstract system (1.1), making precise the conditions on each of the terms and 
reviewing the notions of both martingale and pathwise solutions. We also recall various results from abstract 
probability theory and functional analysis which play a fundamental role in the analysis. 

2.1 Abstract Spaces and Operators 

We begin by fixing a pair of separable Hilbert spaces H Z) V, and assume that the embedding is dense and 
compact. We may thus define the Gelfand inclusions V C H C V, where V' is the dual of V, relative to H. 
We denote by (•, •), | • |, ((•, •)) ^^'^ II ' il the norms and inner products of H, and V respectively. The duality 
product between V' and V is written (•, •). 

2.1.1 The Principal Linear Operator 

Wc now give the precise assumptions on each of the terms appearing in (1.1). They are of course designed 
to include the case of the primitive equations of the ocean, (6.1)- (6.5) as we explain below in Section 6. We 
begin with the linear term supposing that A : D(A) C i? — > -ff is an unbounded, densely defined, bijective. 
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operator such that {AU, [/") = ((f7, [/")) for ah U, [/' E D{A). As such we see that A is symmetric and may 
be understood as bounded operator from V to V' with the duahty product given by 

{AU, U^) = (([/, t/")), for all U eV. (2.1) 

We further sec that A^"^ is continuous as a map from H into V . Since by assumption V CC H , it follows that 
A~^ is compact on H . We may also deduce from the given assumptions on A that A~^ is symmetric. Applying 
the classical theory for symmetric compact operators we infer the existence of a complete orthonormal basis 
{^fe}fc>o for H of eigenfunctions of A so that the associated sequence of eigenvalues {Afe}fe>o form an 
increasing unbounded sequence. For the Galerkin scheme below we introduce the finite dimensional spaces 

Hn = span{$i,...,$„} 

and let Pn, Qn = I ^ Pn be the projection operators onto _ff„ and its orthogonal complement. 

Using the basis {^k} we may also define the fractional powers of A which are also relevant to the analysis. 
Given a > 0, take 

D{A-)^lueH:Y,Xl'^\Uk\'<^\ 
where Uk ~ {U, ^k)- On this set we may define A" according to 

A"U = J2 KUk<^k, iorU = J2 ^i''^''- 

k k 

Accordingly we equip D(A") with the Hilbertian norm 



\UU:=\A-U\=(y^\r\Uk\' 



Classically we have the following generalized Poincare and inverse Poincare inequalities: 

\PnUU, < X'^^-^^PnUU, \QnUU < -^IQ«C/|a., (2.2) 

An 

valid for any ai < a2- 

Note that as in [46] one may verify that D{A^) C D{A") is a compact embedding whenever (3 > a. 
Using (2.1), one may readily verify that D{A^^'^) — V and that \\U\\ = |C/|i/2 for all U €V. Thus, it is clear 
that, in particular, the embedding D{A) C F is compact. 

2.1.2 The Nonlinear Terms 

We turn next to B which we assume to be a bilinear form mapping V x D{A) continuously to V' and 
D{A) X D{A) continuously to H. Furthermore we assume the following properties for B: 

{B{U,U^),U^)=0 ioTa\\UeV,U^ eD{A), (2.3) 

\{B{U,U*),U^)\<co\\U\\\AU^\\U^\\ for all U,U^ eV,U^ eD {A), (2.4) 

\{B{U,U^),U^)\ < co\\Uf^^\AU\'/^\\U^\^^^\AU^\^/^\U^\ for all U,U^ G D{A),U^ e H. (2.5) 

Note that, for brevity of notation, we will sometimes write B{U) for B{U, U). 

We next describe the conditions imposed for F and <t. To this end we introduce some further notations. 
Given any pair of Banach spaces X and y we denote by Bndu{X , y), the collection of all continuous mappings 

* : [0,oo) X Af ^3^ 
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so that 

|l*(x,i)b<c(i + |lx|U), xex,t>o 

where the numerical constant c may be chosen independently of t. If, in addition, 

||*(a;,t) - *(2/,i)b < c\\x~y\\x, x,yeX,t>0 

we say ^ is in Lipu{X, y). 
For F we assume that 

F:[0,oo)xV^H. (2.6) 

In Section 4 we assume that 

FeBndu{V,H). (2.7) 

Further on in Section 5 

FeLipuiV,H). (2.8) 

Similar conditions are also imposed on a. We shall assume throughout this work that 

cr: [0,oo) xi7^L2(H,i/). (2.9) 

Here il is an auxiliary Hilbert space and L2(il, H) is the collection of Hilbert-Schmidt operators between ii 
and H . See Section 2.2 for further remarks. For the case of martingale solutions considered in Section 4, we 
assume that 

a e BnduiH, L2iii, H)) n BnduiV, L2ili,V)) n BnduiD(A), L2iii, D(A))). (2.10) 

On the other hand for pathwise solutions, Section 5, we posit 

a e Ltp,,iH,L2iii,H))nUp^iV,L2i!d,V))nLip^{DiA),L2iii,DiA))). (2.11) 

2.2 The Stochastic Framework 

In order to define the remaining terms in (1.1), that is a(U)dW we must recall some basic notions and 
notations from stochastic analysis. For more theoretical background on the general theory of stochastic 
evolution systems we mention the classical book [11] or the more recent treatment in [43]. 

To begin we fix a stochastic basis S := {il,J^,{J-t}t>a,^,{W^}k>i), that is a filtered probability space 
with {W^'^}/c>i a sequence of independent standard 1-d Brownian motions relative to J^t- In order to avoid 
unnecessary complications below we may assume that J-'t is complete and right continuous (see [11]). Fix 
a separable Hilbert space it with an associated orthonormal basis {efe}fe>i. We may formally define W by 
taking W = J^k ^k&k- As such M^ is a 'cylindrical Brownian' motion evolving over il. 

We next recall some basic definitions and properties of spaces of Hilbert-Schmidt operators. For this 
purpose we suppose that X and X are any pair of separable Hilbert spaces with the associated norms and 
inner products given by | • |x, | • \x ^^'^ {'■• ')x (•, ■)x, respectively. We denote by 



L2iii,X) ^^Re Ci!d,X) : }_^ \Rek\Jc < «^ 



the collection of Hilbert-Schmidt operators from il to X. By endowing this collection with the inner product 
(i?, S')i2(u,x) == J2k(^^k,Sek)x, we may consider L2{ii,X) as itself being a Hilbert space. Note that when 
R € i2(il, X) we shall often denote Rk = Rtk and we may therefore associate R with the sequence {Rk}k>i- 
One may readily show that if i?(i) € i2(il,X) and i?(2) g L{X,X) then indeed i?(2)i?(i) g i2(il,^). 
We also define the auxilary space ito 3 il via 

ilo := < w = ^ akCk ■■ ^ ttfcfc^ 

I fc>0 k 
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endowed with the norm |f |^( :— J2k W^ "^ ~ Sfc ^k^-k- Note that the embedding of it C ilo is Hilbert- 
Schmidt. Moreover, using standard Martingale arguments with the fact that each Wk is almost surely 
continuous (see [11]) we have that, for almost every w € fi, W{oj) € C([0,T],ilo). 

Given an X-valued predictable^ process G G L?{VL]L'^^^{\fi,oo),L2{^,X))) one may define the (Ito) 
stochastic integral 

Mt := / GdW = Y,I GkdWk, 
Jo J. JO 

as an element in A^^-, that is the space of all X- valued square intcgrable martingales (see [43, Section 2.2, 
2.3]). As such {Mt}t>o has many desirable properties. Most notably for the analysis here, the Burkholdcr- 
Davis-Gundy inequality holds which in the present context takes the form. 



E sup 

\ie[0,T] 



GdW 






<cE( / \G\l^iix,x)dt] , (2.12) 



valid for any r > 1. Here c is an absolute constant depending only on r. We shall also make use of a variation 
of this inequality, established in [17] which applies to fractional derivatives of M*. For p >2 and a G [0, 1/2) 
we have 



E 



GdW 





W^-p([0,T]-X}. 




<cE| / lGli,(i^,;,)dt , (2.13) 



which holds for all X-valued predictable G € LP{^;L^^^{[0,oo),L2{ii,X))). For the convenience of the 
reader, we shall recall the definition of the spaces W"'P{[0,T],X) in Section 2.4 below. 

Remark 2.1. Under the assumptions, (2.10), (2.11), on a, the stochastic integral t t^ L <7{U)dW may be 
shown to be well defined (in the ltd sense), taking values in H whenever U G L'^{Q,, Lf^^([0,oo)] H)) and 
is predictable. Such terms may be seen to cover a wide class of examples, including but not limited to the 
classical cases of additive and linear multiplicative noise, projections of the solution in any direction, and 
directional forcings of Lipschitz functionals of the solution. See e.g. [23] for further details. 

In Section 8 we establish the following convergence theorem for stochastic integrals. This result will be 
used below to facilitate the passage to the limit in the Galcrkin scheme. The statement and proof generalizes 
ideas found in [1]. 

Lemma 2.1. Let (51, J^, P) be a fixed probability space, X a separable Hilbert space. Consider a sequence of 
stochastic bases 5„ = (57, J^, {J-'"}t>o,P, VF"), that is a sequence so that each W"" is cylindrical Brownian 
motion (over ilj with respect to J^". Assume that {G"}„>i are a collection of X-valued J^" predictable 
processes such that G" e L'^([0,T], L2{ii, X)) a.s. Finally consider S = (fl, J^, {J^t}t>o, P, VK) and G G 
L'^{[0,T], L2{ii, X)), which is Tt predictable. If 

W^ ^W m probability m G([0, T],ilo), (2.14a) 

G'^^G in probability in L'^{[0,T];L2iii,X)), (2.14b) 

then 

I G'^dW ^ I GdW in probability in L^{[Q,T];X). (2.15) 

JO JO 



^For a given stochastic basis <S, let <t> = f2 x [0, oo) and take Q to be the cr-algebra generated by sets of the form 

(s, i] X F, < s < i < oo, F e JP's; {0} X F, F G J^o- 

Recall that an X valued process U is called predictable (with respect to the stochastic basis S) if it is measurable from (<J>, Q) 
into (X, B{X)), B(X) being the family of Borel sets of X. 
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Finally we describe the assumptions for the initial condition Uq which may be random in general. In 
Section 4, where we consider the case of Martingale solutions, since the stochastic basis is an unknown of 
the problem we are only able to specify Uq as an initial probability measure fiQ on V such that: 

UW^d^iaiU) <oo. (2.16) 

V 

Here q > 2 will be specified below, see Theorem 2.1 as well as Lemma 3.1. On the other hand for pathwise 
solutions where the stochastic basis S is fixed we assume that relative to this basis Uo is 'dV valued random 
variable such that 

Uo e i^(51; V) and is J^q measurable. (2.17) 

2.3 Definition of Solutions 

We next give the definitions of local and global solutions of (1.1) for both Martingale and Pathwise Solutions. 

Definition 2.1 (Local and Global Martingale Solutions). Suppose fio is probability measure on V satisfying 
(2.16) with q > 8 and assume that (2.7) and (2.10) hold for F and a respectively. 

(i) A triple {S,U,t) is a local Martingale solution if S = (il, J^, {J-'t}t>o,P, W^) is a stochastic basis, r is 
stopping time relative to J-t and [/(•) — U{- A r) : fi x [0, oo) -^ V is an Tt adapted process such that: 



U{-At) eL^in;Ci[0,oo);V)) 
Ult<r e L^{n;Ll^i[0,cx,);D{A))) 



2.0..2 (2-18) 



the law ofU{0) is /iq i.e. l^o{E) = ¥{U{0) e E), for all Borel subsets E ofV, and U satisfies for every 
t> 0, 

nt/\T ptAr 

U{tAT)+ {AU + B{U)+F{U))ds^U{0)+ a{U)dW, (2.19) 

Jo Jo 

with the equality understood in H . 

(ii) We say that the (Martingale) solution {S,U,t) is global if t = cxj, ^-a.s. 

We next define pathwise solutions of (1.1) where the stochastic basis is fixed in advance. 

Definition 2.2 (Local, Maximal and Global Pathwise Solutions). Let S = (J7, J^, {Jt}t>o,P, VF) be a fixed 
stochastic basis and suppose that Uo is an V valued random variable (relative to S) satisfying (2.17). As.sume 
that F satisfies (2.8) and that (2.11) holds for a. 

(i) A pair {U,t) is a a local pathwise solution of (1.1) if t is a strictly positive stopping time and U{- At) 
is an Tt-adapted process in V so that (relative to the fixed basis S) (2.18), (2.19) hold. 

(ii) Pathwise solutions of (1.1) are said to be (pathwise) unique up to a stopping time t > if given any 
pair of pathwise solutions {U^ , r) and (C/^, r) which coincide at t ~ on a subset f2 offl, O = {[/^(O) = 
t/2(0)}, then 

P(ljj(C/i(tAr)-t/2(iAr)) =0;Vt>0) =1. 

(Hi) Suppose that {t„}„>i is a strictly increasing sequence of stopping times converging to a (possibly 
infinite) stopping time ^ and assume that U is a predictable process in H. We say that the triple 
{U,^, {■'■ti}ti>i) is a maximal strong solution if ([/, t„) is a local strong solution for each n and 

sup \\U\\^+ I \AU\^ds ^oo (2.20) 

te[o,?] Jq 
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almost surely on the set {^ < 00} . //, moreover 

sup !|C/|p+ / \AU\^ds^n, (2.21) 

«e[o,r„] Jo 

for almost every oj G {^ < 00} then the sequence t„ is said to announce any finite time blow up. 

(Hi) If ([/, C) *'' o, maximal strong solution and ^ = 00 a.s. then we say that the solution is global. 

We may now state precisely the main results in the work: 

Theorem 2.1. 

(i) Suppose that /io satisfies (2.16), for q >S and that F and a maintain (2.7), (2.10) respectively. Then 
there exists a local Martingale solution (5, [/, r) of (1.1). 

(ii) Assume that, relatively to a fixed stochastic basis S, Uq satisfies (2.17) and that F and a fulfill (2.8) 
and (2.11). Then there exists a unique, maximal pathwise solution, (C/, C, {t„}„>i), of (1.1). 

The compactness arguments leading to Theorem 2.1 are carried out in Sections 4 and 5 for (i) and (ii) 
respectively. We provide the details of the passage to the limit needed for both items in Section 7. 

Remark 2.2. (i) We note that, as we are working at the intersection of two fields, the terminology may 
cause some confusion. In the literature for stochastic differential equations the term "weak solution" 
is sometimes used synonymously with the term "martingale solution' while the designation "strong 
solution' may be used for a "pathwise solution". See the introductory text of 0ksendal [38] for example. 
The former terminologies are avoided here because it is confusing in the context of partial differential 
equations. Indeed, from the PDF point of view, strong solutions are solutions which are uniformly 
bounded in H^ , while weak solutions are those which are merely bounded in L^ . In this work we are 
therefore considering both weak and strong solutions in probabilistic sense. From the PDF point of view 
we may say that we are considering strong type solutions since, in the applications considered here V 
is taken to be an appropriate subspace of H^ that incorporates the boundary conditions, etc., for (6.1). 

(ii) The notion of global existence, both for the Martingale and the Pathwise contexts, are included here for 
the sake of completeness. Of course, the passage from the maximal to the global existence of pathwise 
solutions is a significant further step in the analysis and requires further structure for (1.1). We refer 
the reader to [23] and [24] where this is done for the 2D Navier-Stokes Fquations and the 2D Primitive 
Fquations respectively. Current work, making use of the main result herein, treats the global existence 
of solutions for the stochastic 3D Primitive Fquations [13]. 

(Hi) In Section 5, 7 we consider both Martingale and Pathwise solutions of the modified system 

dU + {AU + e{\\U -U4)B{U) + F{U))dt^a{U)dW, U{0) = Uo, (2.22) 

where 

jU, + AU, = Q, U{0) = Uo, (2.23) 

and 9 is a smooth cut-off function as defined below in (3.1). The notions of solutions for (2.22) are, 
with trivial modifications, identical to Definitions 2.1, 2.2 given for (1.1) above. 

2.4 Compact Embedding Theorems 

We shall make use of two compact embedding results taken from [17] which we restate here. See also related 
results in [46]. To this end we first recall some spaces of fractional (in time) derivative. Such spaces are 
natural since we do not expect solutions of stochastic evolution systems to be difFcrentiablc in time but 
merely Holder continuous of order strictly less than 1/2. 
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Let X be a separtable Hilbert space and denote the associated norm by | • \x- For fixed p > I and 
a G (0,1) we define 



W'^^P{[0,T];X) := lu G L^{%nX), j^ j^ ^E^l^J^^^Ldt'dt" < ^ 



f^ /-^ \TT(i'\ _ TT(+'I\\P 

/o Jq 
We endow this space with the norm 



For the case when a = 1 we take W^'P{[0,T];X) := {U G LP{[0,T];X);^ G LP{[0,T];X)}, to be the 
classical Sobolev space with its usual norm 



\U\w^.P{[O.T]:X) '■— / |f^(OI 



/ ""v + |f<''> 



p 



Note that for a G (0,1), W^'P{[0,T];X) c VF"-P([0,i];X) and |f/|w'°.p([o,T];X) < C'|?7|,^i.p([o,t];X)- With 
these preliminaries in hand we may now state the compact embeddings needed below (see [17]) 

Lemma 2.2. 

(i) Suppose that X2 Z) Xq D Xi are Banach spaces with X2 and Xi reflexive, and the embedding of Xi 
into Xq compact. Then for any 1 < p < 00 and < a < 1, the embedding: 

LP([0,r];Xi)nW^"^P([0,r];X2)ccLP([0,r];Xo) (2.24) 

is compact. 

(ii) Suppose that Yq 1) Y are Banach spaces with Y compactly embedded in Yq- Let a G (0, 1] and p G (1, 00) 
be such that ap > 1 then 

W°''P{[0,T];Y) CC C([0,r],ro) (2.25) 

and the embedding is compact. 

2.5 Some Tools Prom Abstract Probability Theory 

We next review some classical convergence results for probability measures defined on separable metric 
spaces. In conjuction with the embeddings given in Section 2.4, these results provide some powerful means 
to address the difficulty of establishing compactness for the collection of Galerkin approximations associated 
to (1.1). 

Let {X,d) be a complete separable metric space and take B{X) to be the associated borcl cr-algebra. 
Also, we define Cb{X) to be the collection of all real valued continuous bounded functions on X and take 
Pr{X) to be the set of all probability measures on (X, B{X)). Recall that a collection A C Pr{X) is said to 
be tight if, for every e > 0, there exists a compact set K^ (Z X such that: 

M(-^e) > 1 - e for all /i G A. 
On the other hand a sequence {/in}ri>o C Pr{X) is said to converge weakly to an probability measure /i if 

fd^in -^ I fdfj. 



over all / G Cb{X). We say that a set A C Pr{X) is weakly compact if every sequence {p,n] C A possesses 
a weakly convergent subsequence. 

Proofs of the following classical results may be found in e.g. [11]. 
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Proposition 2.1. (i) A collection A C Pr[X) is weakly compact if and only if it is tight. 

(ii) Suppose that a sequence {^n\n>i converges weakly to a measure fi. Then there exists a probability 
space {Q,T,¥) and a sequence of X valued random variables {Yn}n>o (relative to this space) such that 
Yn converges almost surely to the random variable Y and such that the laws of Yn and Y are fin and 
II, i.e. fin{E) = P(F„ e E), fi{E) = ¥{Y G E), for all E e B{X). 

Finally wc come to an elementary but powerful characterization of convergence in probability introduced 
in [28]. Suppose that {F„}„>o is a sequence of X-valued random variables on a probability space (17, J^, P). 
Let {/XTi.rn}ri,Tn>i be tlic coUcction of joint laws of {Yn\n>i, that is 

/i„,„(S) := P((y„,y™) e ^), £ e B(X X X). 

The result from [28] is the following: 

Proposition 2.2. A sequence of X valued random variables {Yn\n>o converges in probability if and only 
if for every subsequence of joint probabilities laws, {/infc,mfc}fc>0; there exists a further subsequence which 
converges weakly to a probability measure fx such that 

n{{{x,y)eX X X -.x^y})^!. (2.26) 

3 The Approximation Scheme 

We now implement a Galerkin scheme for (1.1). To this end we introduce the projected operators 

S"((7) = PnB{U), F^{U) = PnF{U), a^{U) = P„a([/), 
where U ^ V. We shall also make use of a 'cut-off' function : M ^- [0, 1], which is C°° and such that: 

, , I 1 if jx'j < K, , , 

Here we choose k to be any positive constant, independent of n, such that 

^<^, (3.2) 

64co 

where cq is the constant appearing in (2.5). The reason for this choice will be made apparent in the proof of 
Lemma 3.1 (see (3.19), (3.26)). 

Wc now fix a stochastic basis S = {fl,J-,{J-t}t>o,^,W) and an element Uq ^ V with law fiQ. We find 
pathwise solutions to the Galerkin systems defined by (3.6) relative to this basis below. Since we allow for 
an ill-behaved nonlinear term B (that is satisfying (2.3), (2.4), (2.5)) wc introduce an auxiliary linear system 
in order to carry uniform estimates on the Galerkin systems. We take C/" to be the unique (i/„ valued) 
solution of 

^u: + Au:: = o, kio) = p„[/o. (3.3) 

at 
One may readily verify that, for any p> 2, [/" satisfies the estimates 

sup \\uiT + r \Au:mu:r-'dt' +{ I \Au:\^dt'\ <c\\uor. (3-4) 

t'<t Jo 




It is also clear that 

\U:\w^.^o,T,H) < c I \AU:\'dt' < cWUof. (3.5) 
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With these notations in place we define the Galerkin system at order n 

(3.6) 



dC/" + [AC/" + 0(||t/" - C/,"||)B"([/") + F'^{U'^)]dt = cr"(C/")dW^, 
C/"(0) = P„C/o := (7o". 



Here C/" is an adapted process in C([0,T];i7„) = C([0,T],R"). The C/" appearing in the cutoff function 
are solutions of the linear systems (3.3). The significance of this addition will become clear in the proof 
of Lemma 3.1 below. Note that, due to the preserved cancellation property in the nonlinear portion of the 
equation, the existence and uniqueness of solutions at each order is standard. See, for example, [16] for 
further details. 

Lemma 3.1. Assume that F and a satifisfy (2.7), (2.10). Fix Uq, a V-valued, Fq measurable random 
variable and consider the associated sequence of solutions {C/"}n>i of the Galerkin system (3.6), (3.3). We 
suppose that the constant k appearing in the cutoff function 9 satisfies (3.2). Let p > 2 and suppose that 



t/o||' < oo for some q > max{2j5, 4}. 



(3.7) 



Then there exists a finite number K > depending only onp, Ej|C/o||'' o,'>^d the the constants in (2.5), (2.7), 
(2.10) such that 



(i) for every n > 1, 



and also 



and finally 



E sup |lJ7"f + / |AC/"n|C/"||P-2^i' < K, 

\t'<T Jo / 






|Ac/'Tdt' 







E 



cr"(C/")dl^ 



< K. 



(3.8) 



(3.9) 



(3.10) 



(ii) If under the given assumptions we additionally suppose that p > 4, then we have, for all n>l: 

2 \ 



U"{t)- I cr''{U")dW 




<K. 



(3.11) 



Wi-2([0,T];H)^ 



Proof. Define C/" := t/" — t/". We may readily observe that t/" satisfies 

dtJ" + [AC7" + 6i(||C7"||)B"(C7" + U^) + i^"(t7" + U^)]dt = a"{U" + U^)dW, 

C7"(0)=0. 

We apply A^^^ to this system. With the Ito formula we infer, for p > 2 that, 

d||c7"f +p|Ai7"H|i7"f-2dt 

+ ^^^^^—^{a'^iU" + U:), AU'')^\\U''\\P-^dt - pe{\\U"\\){B'^U^^ 

+p||f/"f-2(o-"(i7" + [/;), Ai7")dp^ 

=:( Jf + Jf + JP + Jl)dt + jPdW. 



(3.12) 



(3.13) 
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We are able to estimate the first four deterministic terms pointwise in time. Using (2.7) we observe that 

I j[i <c(i + \\u:\\ + iic7"ii)|Ac7"iiif7"r-2 < i\Au-\^u"\r' + c(i + \m\\ + iic/"ii)'iic/"r-- 



<f I^c7"n|i7"f -2 + c(i + \\U'^\\)P + c||i7"||P. 

8 
The terms J2 and J^ are also estimated directly using (2.10) 

kf I + \Ji\ <c((i + mi)2 + \\u"r)\mr' < c(i + \\u:\\r + c^-r. 

Using the bilinearity of B the term Jf splits according to: 

I j|i <p^(iit7"ii)i!^"r-'iA;7"i(|i?([/:)i + \B{u", u:)\ + \b{u:, u")\ + \b{u^)\) 

• — '^4,1 + ^^4,2 + ^^4,3 + '^4,4- 

We estimate each of these terms using (2.5). For J^^ we have 

' 8 

<^\\U"\\P~^\AU"\^ + cKP-^\\Ul'f\AU:\^ 
8 

For the next two terms we estimate 



(3.14) 

(3.15) 
(3.16) 

(3.17) 



^4,3 1 



<^||^"f-2|AC7"|2 + c6l(||t/"||)||C7"f ||C/;||2|yl[/:'|2 
8 

<^||^"f-2|^i7n|2^_c^P||[/n||2|^^n|2^ 

8 
The last term yields to the bounds 

1 44I < cop6'(||C/"||)||L/"f-VC^"P < 2kpco\\U"\\p-^\AU"\^ < ^\\U"\\p-^\AU''\'^. 
' 8 

Note here that the last inequality follows from the requirement (3.2) imposed on k. 

Finally we address the stochastic terms. Observe that for any pair of stopping times < Tq < rt < T, 
the BDG inequality, (2.12) with r = 1, implies that 



(3.18) 



(3.19) 



E sup 

Ta,<t<Tt 



J^dW 



< cE 



< cE 



' \\U'''\\^^P-^\a'''{U'' + U^),AU")^d! 



1/2 



||^n||2(p-l)(^^||[;„||_^||j^n||^2^^ 

<c(e sup \\u''r-'( r\i + \\u:\\ + \\u^)^ds) 



1/2^ 



(3.20) 



<1e 

- 2 

<iE 

- 2 



sup iic7"r) +cE f rv + \\u:\\ + \m\)'ds^ ' 

Ta<t<Tb / \Jr„ 

sup ||C7"r) + cE /"'((I + \\U:\\)P + \\U'T)di 

.Ta<t<Tb I Jva 



Combining the estimates (3.14), (3.15), (3.16), (3.17), (3.18), (3.19), (3.20) we find, for any t e (0,T], 

¥.( sup ||t/"r+ / |^C7'f||t/"f-2dt'^ <cE /" (||C/"r + (l + |ylC/:f ||C/r|p + ||C/rir))di' 
\t'e[o,t] Jo / Jo 



<c 



: sup iic7"r+E(i + |Ac/,"njt/,"f + iit/:r) ut'. 

se[o,t'] / 



(3.21) 
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Applying then the Gronwall inequahty yields 

e( sup ||C/"|r+ / |^C7"|2||c7"f-2rfi' I <cE / (l + |^C/:f ||C/r|p + ||C/rir)di' 



,0<t'<T 











(3.22) 



<cE(l + i|[/o||)'"^^<P'''^ 



The second inequahty follows from (3.4). We also note that the term involving |^C/"p||[/"|p is responsible 
for the first part of the moment condition (3.7). 

In order to complete the proof of (3.8) we observe that 



7. \ p/2N 

n|2 ../ 



E( sup ||C/"f + / \AU"\^\\U"\\P-'^dA <cE\ sup ||C/"f+l / |ylt/"|^df 

Vo<t'<T Jo J \o<t'<T VJo 



<cE I sup ||C/:'f - 

0<t'<T 






\AU':\^d£ 



cE I sup ||C/"||P- 

0<t'<T 






(3.23) 







Given the estimates (3.22) for [7" and (3.4) for L/" it therefore remains to estimate the last term, i.e. to 
prove the analogue of (3.9) for C7. Returning to (3.13) for the case p = 2 we must therefore find suitable 
estimates for the left hand side of the expression 

p/2 



E / \A\JYdt 



p/2 



<E 



( / \Jl\ + \Jl\ + \4\ds+ sup 
\Jo te[o,T; 



JldW 



(3.24) 



/o / \Jo telo.T] 

Note that, when p = 2, | J|| = 0. By treating \J'f\, | J|| in a similar manner to (3.14), (3.15), we infer 

I Jl I + I J|| < 2-(2+2/p)|^i/n|2 ^ ^(1 ^ j|[;nj|2) ^ c!|C7"||2. (3.25) 

For I Jf I we estimate similarly to (3.16), (3.17), (3.18), (3.19) to deduce 

|j2| <2-('^+^/p)\AU''\'' +c\\u:\\^\AUl'\^ +Akco\AU\'' <2-^^+^/P^\AU''\^ + c\\u::f\AU^\''. (3.26) 

The constant cq after the first inequality is from (2.5). Thus, the assumption (3.2) justifies the second 
inequality. For the stochastic intergral term in (3.24) we apply the BDG inequality, (2.12), and deduce: 



E sup 
te[o.T] 



JldW 







p/2 



<cE ( j {a''iU'' + Ul'), AU")'^dt' j < cE 

<c^(J {^ + \\u''\\^ + \\u:T)dA <ceJ (i + iic/"r + iic/,"r)di' 






(3.27) 



Applying (3.25), (3.26) and (3.27) to (3.24) we have 



f-T _ \ p/2 

E( / \AU\^dt' 



<-E 



p/2 



\Au\'ds\ +cE/ (i + i|[/:r + ii(7"r)di'+cE / \\u:\f\Au:fds 



p/2 



(3.28) 



p/2 



<-E / \AU\^ds] +cE sup (l + ||C/:r + |lC7"r)+El|L/ofP. 



te[o,T] 
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Note that the terms involving ||f7"|p|y4t/"p are treated in the final inequality using (3.4) and are responsible 
for the second part of the moment condition (3.7). Applying (3.28) in turn to (3.23) we finally conclude 
(3.8). With (3.4), (3.9) also now follows from (3.28) 

The bound (3.10) is a direct application of (2.13) with (2.10): 



CT"(C/")dM^ 



<cE / l-"^rr»MP 
Jo 



W''iUn\UuM)dt<cE il + \UY)dt. (3.29) 



We finally establish (3.11). Integrating (3.6) we observe that 



[/"(<) - / a"([/")dl^ ^K+ I [AV^ + 6'(|1C/" - C/,"||)B"(C/") + F"(;7")]rf<. (3.30) 

Jo Jo 



2 7^ 

<c\Uo? + c / (|^C/"|2 + |B"(f/")|2 + |F"(t/")|2)ds 

W^-'^{[0,T]-H) Jo 

<c\Uo? + cl (|^C/"|2 + |B(C/")p + |F(C/")nds (3.31) 

Jo 

1 2 I ^ / /I I l|rrM||2\/-i I I ^7-7-n|2\ 



With, (2.7), (2.5) wc infer: 



f/"(t)- / fj"(f/")dT^ 





<c\UoV + c (1 + l|C/"lr)(l + \AU''\')ds. 
Jo 

Taking expected values in this expression and applying (3.8) (i) for the case p = A gives (3.11). The proof 
is now complete. D 

4 Local Existence of Martingale Solutions 

In this section we establish the existence of a Martingale solution of (1.1). The first step is to make use of the 
uniform estimates established in Lemma 3.1 we infer the compactness (in certain spaces) of the probability 
laws associated to the Galerkin approximations. We then change the underlying probabilistic basis in order 
to find a new sequence of random elements equal in law to the original Galerkin approximations but which 
converge almost surely. The technical details of the passage to the limit, which is used also below for the 
case of pathwise solutions, is carried out in Section 7 below. 

4.1 Compactness Arguments 

For a given initial distribution /ig on V wc fix a stochastic basis S = {^,J-^ {J^t}t>o,P, W^) upon which 
is defined an J-q measurable random element t/g with distribution /ig- Consider the sequence of Galerkin 
approximations {[/"} solving (3.6) relative to this basis and initial condition. Wc consider the phase space: 

A't/ = (L2(0,T;^)nC([0,r],\/'), A'M/ = C([0,T],ilo), A" = A-y x A-^. (4.1) 

We may think of the first component, Xij, of this phase space as the set where the solution [/" lives and the 
second component, Xw^ as being the set on which the driving Brownian motions are defined. We consider 
the probability measures 

M&(-) = P(C/" e •) e Pr(L2([o,T];y)nC([o,r],y')), (4.2) 

and 

m(-) = ^^w{■) = nW e •) e Pr(C([0,T],ilo)). (4.3) 

This defines a sequence of probability measures /i" — fi'jj x /^^^ on the phase space X. By making appropiate 
use of Lemma 3.1 we will now show that this sequence is tight. More precisely: 
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Lemma 4.1. Suppose that /ig satisfies (2.16) with q > 8. Consider the measures fi" on X defined according 
to (4.2), (4.3). Then the sequence {/i"},i>i is tight and therefore weakly compact over the phase space X . 

Proof. By applying Lemma 2.2, (i) with X^i ~ H, Xq ~ V, Xi ^ D{A), p = 2 and a = 1/4 we deduce that 

L^{[0,T];D{A))r\W^/^'^i[0,T];H)ci(lL^{[0,T];V). 

For i? > we define the set 

B},^{Ue L\[0,T];D{A)) n W'/'^'{[0,T];H) : |t/|i.([o,T];D(A)) + \U\W.HIo,t];H) < ^'} 

which is thus compact in L'^{[0,T],V). Due to the Chebyshev inequahty and the uniform estimates (3.8), 
(3.11), (3.10) in the case p = 2, we estimate, 



f^uii^R) ) — IP'(I^"Il2([0,T]:D(A)) + l^"lwi/4,2([o^T]:H) - ^ ) 

c 



(4.4) 



2 

-i?2 



lAV^l^dt' + \U''' 



'wt'^{[o,t]-h} 



< 



i?2' 



where the numerical constant c is independent of n. 

Choose a € (1/g, 1/2) so that aq > 1. By Lemma 2.2, (ii) with Ya ^ V ^ D^A'^^'^) and Y = H we 
infer the compact cmbcddings 

W^^^{[0,T];H) CC C([0,T],F'), W°'''^{[0,T];H) CC C([0,T], F')- 

For i? > 0, we take B|'^ and B"^/ to be the balls of radius R in W'^'^{[0,T], H) and VF"''?([0,r], iJ) 
respectively. It follows that for R> 0, Bj^ := Bj^ + Bj^ is compact in C([0,T], V). Since indeed, 

{[/" G B|} D |t/"(i) - / a''{U'')dW e B^;H n | / cr"([/")dw^ e B^/\ , 

we may apply Chebyshev's inequality and then the uniform estimates (3.11) (3.10) to infer 

2 \ 



/x^((B^)'-)<P 



t/"(i)- / fT"(C/")dW^ 





> i?^ 



fT"(C/")dM^ 



vyi.2([o,T];_ff) 



(4.5) 



VK°''J([0,T];i/) 



<- 



■R^' 



As above the c is independent of n. 

It is not hard to see^ that B}^ n Bj^ is compact in L'^{0,T;V) n C{[0,T],V') for every i? > 0. As a 
consequence of (4.4) and (4.5) we have 



f^uUB'n n Bj,f) < fi-uUB'nf) + l^l{{Blf) < 



i?2 



We therefore take A^ := B^ . ni?^ . , with c the constant which appears on the left hand side immediately 

-v/2c/e -v/2c/e 

above. With this definition we infer that for e > 0, 



y^UAe) > 1 



(4.6) 



20nc need only verify that if {(7"}„>o C ^^(o, T; V) n C([0, T], y' ) and if U" ^ U in L2(o, T; V), (7" -> (7 in C{[0, T], F') 
that U = U 
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over all n. 

We next turn to the sequence {iJ,^}n>o- This sequence is constantly equal to one element and is thus 
weakly compact. Hence, as a consequence of Proposition 2.1, (i) {/x^}n>o must be tight. We therefore infer 
the existence of collection of compact sets A^ c C([0,T],ilo) so that 

MlV(^e)>l-| (4.7) 

for all n. 

We now have everything in hand to conclude the tightness and therefore the weak compactness of {/x"}„>o- 
For e > we define /Cg := A^ x A^ which are compact in X. By (4.6) and (4.7) we infer that, for any e > 
and every n, 

and thus that {/^"}„>o is tight in X . Prohorov's theorem, given herein as Proposition 2.1 therefore implies 
that /J," is weakly compact. The proof is therefore complete. D 

4.1.1 Strong Convergence on the Skorohod Space 

Given ^o (satisfying (2.16) with (j > 8) we have shown that the sequence of measures {/i"'}ra>i associated to 
the Galerkin sequence (C/", W) are weakly compact on X . Passing to a weakly convergent subsequence /x"*" 
we now apply the Skorohod embedding theorem, Proposition 2.1, to infer the following Proposition. 

Proposition 4.1. Suppose that fiQ is a probability measure on V satisfying (2.16) with p > A. Then there 
exists a probability space (r2,J^,P) and a subsequence Uk and a sequence of X valued random variables 
(i7"SW"'=) such that 

(i) ([/"*=, VF"*") converges almost surely, in the topology of X , to an element {U,W). 

(a) W"''' is a cylindrical Wiener process, relative to the filtration J-'^'' , given by the completion ofa(W™''{s), 
C7"'=(s);s<i). 

(Hi) Each pair {(/"■'' ,W"'') satisfies 



dij"" + [AU'"' +6'(||i7"'= - U^^^IDB"" {U"") + F'"' {ij">')]dt = a"" {U"'')dW'"' , 



(4.8) 



where we define U"'' by: 



^ u:' + Au:" =0 u:'{o)^ u^" . (4.9) 



With this proposition established the existence of a local Martingale solution follows once we have shown 
that {U,W) and an appropriately defined stopping time t (see (7.5)) satisfy (1.1). This passage to the limit 
argument, which is technical and delicate, is carried out in Section 7 below. 

Remark 4.1. While Proposition 4-1, ft) follows directly from Proposition 2.1, (ii) further steps are required 
to establish (ii), (Hi). These technical points may be demonstrated in a similar manner to previous works. 
See [1]. 

5 Local Pathwise Solutions 

We turn now to study Pathwise solutions of (1.1). Here the key step is to apply Proposition 2.2 in order to 
show that (U^-jW) converges almost surely in L'^{[0,T];V) n C{[0,T],V') relative to the initial stochastic 
basis. The diagonal condition, (2.26) translates to a question of pathwise uniqueness which we address first. 
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5.1 Local Pathwise Uniqueness 

The following proposition establishes the uniqueness, pathwise, for any pair of solutions of the modified 
system (2.22). Such solutions appear in an intermediate step in the compactness arguments in Section 5.2 
below. 

Proposition 5.1. Suppose that (S, t/'^') and (S, t/"') are two global Martingale solutions of (2.22) relative 
to the same stochastic basis S := {n,T,{Tt}t>o, P, VF). Assume that, in addition to the conditions imposed 
in Definition 2.1, F and a satisfy the Lipschitz conditions (2.8) and (2.11). Define 

fio = {C/^'H0) = t/(2)(0)}. (5.1) 

Then U^^' and U^"^' are indistinguishable on f^o in the sense that 

P (loo(C/'(t) - U^{t) = 0;Vt > 0) = 1. (5.2) 

Remark 5.1. We note that, with trivial modifications to the proof that follows, one may establish that 
Pathwise solutions of (1.1) are unique in the sense of Definition 2.2, (ii). 

Proof of Proposition 5.1. Define R ~ U'^^^ — C/^^^ and let R = l^oR- Note that, by definition, R e 
C([0, oo);y) n Lf^^{[0,oo);D{A)), a.s. Due to the bilinear term B, when we attempt to estimate R, stray 
terms arise that involve only [/'^-'^^ or U'^^\ See (5.8), (5.9) below. To remedy this situation we define the 
stopping times 

t(") := inf I r ||C/W||2|A[/(i)|2 + ||(7(2)||2|AC/(2)|2^g > „| (53) 

Clearly this is an increasing sequence. Futhermore, since U^^\ U^'^' are global solutions, we may infer that 
Zz7TT-„_i.oo ''■'•"'' = 00 from (2.18). Hence, the desired result will follow if we show that for any n, T > 0, 

e( sup \\R\\A =0. (5.4) 

\[0,r(")AT] / 

Subtracting the equations (c.f. (2.22)) for U^^^ from that for U'^^^ we arrive at the following equation for 
R: 

dR+iAR + 9{\\U'^'^ - C/i''||)B(L/(i)) - 0{\\U''^'> - C/i'^||)B(t/(2)) 

+ F(C/(i)) - FiU^^^))dt = (a(t/(i)) - a(C/(2)))dM/, (5.5) 

i?(0) = C/(i'(0)-t/(2)(0). 

Ito's lemma yields the following evolution equation for ||i?|p: 

d\\Rf + 2\AR\^ =2(0(|l[/(2) - C/i2)||)B([/(2)) - 0(||C/(i) - ui^^\\)B{U^^^),AR) 

+ 2(^(1/(2)) - F([/(i)), Ai?)di + ||a([/(i)) - <7{U^'^)\\l,^^y)dt (5.6) 

+ 2(cr(C/(i)) - a([/(2)), AR)dW. 

Fix n and stopping times Ta,Ti,, such that < Tq < t;, < t'"'. Integrating in time and taking supremums. 
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multiplying by Iq,, and finally taking an expected value we arrive at the expression 
e( sup ||i?|p+ [ " \AR\^ds] 

<E\\R{Ta)\\^ + 2E r \m\\U^^'> - UP\\) - 0(||(7(2) - (7?^!|))B((7(i)), Ai?)|di 

+ 2E / ''|(B([/(i))-S(C/(2)),Ai?)|df + 2E / " \{F{U'-^^) ~ F{U^^^),AR)\dt ^^■'^^ 



{a{U^'^)-a{U''^'^),AR)dW 



2E sup 

\R{TaW + Jl + J2 + J3 + J4 + J5- 



+ E / InoMU^'^) - <y{U^^^)\\Uny)ds 



For Ji, we use that 9 is Lipschitz. Sec (3.1). Applying (2.5), we have 

Ji <cE / ' ||(L/(i) - [/(')) - (C/i'^ - ui^^)\\\{B{u''^^),AR)\dt < cE /" ' \\R\\\\u'^'^\\\AU'^'^\\AR\dt 

<-E / \AR\^ds + c¥. ||[/(i)f|^C/(i)p||i?||2ds. 

Note that since both t/* , C/* satisfy the linear equation (2.23) it is clear that, for every t >Q lsio(t/* {t) — 

(2) 

[/* {tj) = almost surely. For J2 the bilinearity of B and (2.5) imply: 

J2 = 2E / ' |(S(C/(i) - C/(2), t/(i)) + B{U^^\U^^^ - U^^^),AR)\dt 



<2E / \{B{rM^^) + B{U^'^\R),AR)\dt 

;; (5.9) 

<CE / (||;7W||l/2|AC/(l)|l/2 + ||[/(2)||l/2|^C/(2)|l/2)P||l/2|^^|3/2^^ 
<iE /" V^l'cis + CE /" "(llt/^l^f |A[/(l)|2 + j|[/(2)||2|ylC/(2)|2)p||2^^^ 

The terms J3 and J5 are estimated directly making use of (2.8) to infer 
h <cE I ' loo|i^(C/('') - F([/(2))pi?|ds 

<cE I " \\R\\ \AR\ds < ^E /" ' \AR\^ds + cE f ' \\R\\^ds, 
and making use of (2.11) to deduce 

Jg < cE / \\R\^ds. (5.11 



(5.10) 



Finally, J4 is addressed using (2.12), with r ~\ and then (2.11) 

rn \ 1/2 



J4 <ce( r lo„(c7(L/(i)) - (t([/(2)), Ai?)2dsj 

1 /9 

<cE (^£' T^noMU^'^) <^{U^'^)\\Uu,vm'd^) (5-12) 



Tb _ \ 1/2 1 _ I^Th 

<cE[ I ||i?||M < -E sup ||i?f + cE/ |ji?fds. 

/ 2 tg[r„,^j,] J^^ 
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Applying the estimates in (5.8), (5.9), (5.10), (5.11), (5.12) to (5.7) we infer 



e( sup ||i?|p+ [ " \AR\^ds 

<cE\\R{Ta)f + CE P(|jC/(l)||2|AC/^l)p + ||{/(2)||2|^[/(2)|2_^_i)P||2^^ 



te[r..r,] „,„ . ^^^^g^ 



With this estimate we may finally apply the stochastic Gronwall lemma, as in [23] to conclude (5.4). The 
proof is complete. D 

5.2 Compactness Revisited 

We return to the sequence {t/"} of Galerkin solutions of (3.6) defined relative to the given stochastic basis S. 
We assume throughout this section that Ej] UaW^ < oo for some q> 8. Once we have established the existence 
of local pathwisc solutions for all initial data in this class the general case, (2.17) may be established via a 
localization argument. See e.g. [23]. 

In pursuit of Proposition 2.2 we consider the collection of joint distributions /i^'™ given by ([/",t/™). 
For this purpose we define the extended phase space (cf. (4.1)) 

^ (5.14) 

Xu = L'{^,T-v) n c([o,r],i/'), Xw = c([o,r],Jio). 

As above in (4.2), (4.3) we let \jJ^{E) = P(C/" G E) for E (^ Xu and [iw{E) ^ P(P^ e E) for E e Xw Take 

M^'" = Ma X /i^\ t."-" = M^XM[7xMH'. (5.15) 

Similarly to Lemma 4.1 we prove: 

Lemma 5.1. Suppose E||[/o|l'' < oo for some q > 8. The collection {z/"'"'} (and hence any subsequence 
|jynfc,™fcj.j jg tight (and hence compact) on X^ . 

Proof. The proof is nearly identical to Lemma 4.1. We determine the sets B\(, Bj^ exactly as previously. With 
trivial modifications (see (4.6) and remarks immediately above) we can therefore choose A^, A^ compact in 
A';/ and AV respectively so that /ij}(Ae) > l — j, and iJ.^{A^) > 1—|, for every n. Taking /C^ := A^xA^xA^, 

which is compact in X'^ we see that i'"'"'(/Ce) > (l — ■§) (l — §) > 1 — e, which holds for every < e < 1 
over all m, n. The proof is complete. D 

Suppose now that {h'^'"'^''}k>o is any subsequence. By Lemma 5.1, {i^^'°''"'°}fc>o is tight and hence by 
Proposition 2.1 (i), we may choose as subsequence k' so that v'^k^rn^. converges weakly to an element v' . 
By applying Proposition 2.1, (ii) we next infer the existence of a probability space (17, J^, P) upon which 
are defined a sequence of random elements (t/"'=, C/"*'', W^ ) converging almost surely in X'^ to an element 
{U,ij,W) in such a way that P ('(C/"'s ^7™'^ 1^'=') € ■) = v<^"''>^{-) and fUu.U.W) S •) = v {■). Let 

Zk> = {U<,W^'), Zk' = {U<,W^')., Z ^ {U, W) and Z = {U, W). Note that in particular filj"'""- converges 
weakly to the measure iiu defined by 

^lu{■)■.= niU,U)e■) (5.16) 

Exactly as for Proposition 4.1, we may establish the conditions for Proposition 7.1 below for both Zk/, 

Z and Zk', Z. As such we infer that both U and U are global Martingale solutions of (2.22) over the same 
stochastic basis S = {fl, T , {J^t}t>o, P, W^. Since it may be readily shown from the above convergences that 

f/(0) = t^(0) a.s. we infer from Proposition 5.1 that {/ = [/ in Xjj a.s. In other words 

^{{(x, y)^X^jxX^j:x = y\)= P(L/ = U in Xjj) = 1. 
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With this conclusion, Proposition 2.2, now inipHes that the original sequence C/" defined on the initial 
probability space (fi,^, P) converges to an element U, in the topology of Xjj- By a final application of 
Proposition 7.1^ below we may infer that [/ is a global pathwise solution of (2.22). Hence, taking r as in 
(7.5) below we conclude that {U, r) is a local pathwise solution of (1.1). 

The passage from a local to a maximal pathwise solution in the sense of Definition 2.2, (iii), may now be 
carried out as in [26]. See also [29]. This completes the proof of Theorem 2.1. 

6 Example: The Primitive Equations of the Ocean 

As discussed in the introduction the primary motivation in the development of the abstract theory was to 
be able to treat the existence of local, pathwise solutions of the stochastic primitive equations of the ocean 
considered in the /3-plane approximation. We now recall this system of equations and show how Theorem 2.1 
applies to these equations. 

The stochastic primitive equations of the ocean in their /3-plane approximation take the form 

^tv + (v V)v + wa^v H Vp + /k X V - (UvAv - !/v5^^v = F^ + av{^^,T, S)Wi, (6.1a) 

Po 

dzp=-pg, (6.1b) 

y ■v + dzw = 0, (6.1c) 

dtT + (v • V)T + wd,T - ^itI^T - vrd^^T = Ft + ari^^, T, S)W2, (6. Id) 

dtS + (v • V)^ + wd,S - fisAS - usd,,S = Fs + ^^(v, T, S)W3, (6.1e) 

p = po(l - MT - T,) + PsiS - Sr)). (6.1f) 

Here, w^T, S,p^ p represent the horizontal vclocity,temperature, salt concentration, pressure and density of 
the fluid under consideration; p^, v^. pT, ^t, Ps- ^s arc (possibly anisotropic) coefficients of the eddy 
viscosity and of the heat and salt diffusivity respectively; / is the Coriolis parameter appearing in the 
antisymmetric term in (6.1) and accounts for the earth's rotation in the momentum budget. The stochastic 
terms are driven by white noise processes Wj and are understood in the Ito sense. The equations as given 
above model oceanic flows, however equations of a quite similar structure may be given that describe the 
atmosphere and the coupled oceanic atmospheric system, or the same equations on the sphere. 

Of course (6.1) is supplemented with appropriate boundary conditions which, among other considerations 
must account for the coupling at the oceans surface with the atmosphere. The evolution occurs over a 
cylindrical domain A4 ~ Mq x {—h, 0), where A^o is an open bounded subset of M^ with smooth boundary 
dMf). We denote by uh the outward unit normal to dMo The boundary dA4 is partitioned into the top 
r.j = Mq X {0}, bottom Fb = Mo x {-h} and sides P; = dMo x (-/i,0). Wc prescribe (see [33], [42]) 

i/v^^v + QfvV = 0, w = i^rdzT + arT ^ 0, 0^3 = 0, on P„ (6.2) 

v = 0, w = 0, d,T = 0, d.S^O, on Pft, (6.3) 

v = 0, 5n„r = 0, 5n„5 = onP,. (6.4) 

The equations and boundary conditions (6.1), (6.2), (6.3), (6.4) are supplemented by initial conditions for 
V, T and S, that is 

V = vo, T = To, S = So, at t = 0. (6.5) 



^Actually, in contrast to previous cases above, the convergence is more straightforward since in this case we need to consider 
only one fixed driving brownian motion W throughout. 
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With the boundary conditions, (6.2), (6.3), (6.4), we may reformulate (6.1) according to 



1 r° 

dt^ + (v • V)v + u)(v)5,v + —Vps - g / {Pt'^T + psVS) dz 

Po Jz 



(6.6a) 



+ /k X V - MvAv - i^^dzzw = Fv + av(v, T, S)Wi, 

V • vdz, / V • vdz = 0, (6.6b) 

StT + (v • V)T + w{v)dzT - /irAT - i^t^^^T = Ft + crT(v, T, S')1^2, (6.6c) 

ftS* + (v • V)5' + w{v)dzS ~ fisAS - lysdzzS ^ Fs + as{v, T, S)W3, (6.6d) 

which is the basis for the functional framework that we next recall. Our presentation and notations closely 
follow the recent survey [42] which covers the deterministic setting. 

We will denote by U the triple of prognostic variables, U ~ (v,r, S) (comprising four scalar variables), 
and we set 

H = \{v, T, S) e L^{Mf : V • / vdz = over Ma, nn ■ f vdz ^ over dMo, f SdM = 1, 

I J-h J-h Jm J 

which wc equip with the classical L^ inner product*. Define Ph to be the Leray type projection operator 
from L^(A^)'* onto H. For iJ^(A^)* wc consider the subspace: 

V^ l{v,T,S)eH\M)^ -.^7 ■ f vdz = OovcT Mo,v = on Ti Li Tb, f SdM = o\. 
We equip V with the inner product 

(([/,[/«)) :=((v,v«))i + ((r,r«))2 + ((5,5«))3, 

((v, v))i := / (/^vVv • Vv« + i/v^zv • ^z^r^) dM + a^ [ vv^dT,, 
Jm Jvi 

((T, T»))2 := / {^JiTVT ■ Vr« + VTdzT ■ d^T^) dM + ar [ TTdT^, 
Jm Jvi 

((5, 5«))3 := / {^is^S ■ V5« + ugdzS ■ d^S^) dM. 
Jm 

Note that under these definitions a Poincare type inequality \U\ < C\\U\\ holds for all U gV. We take V(2) 
to be the closure of V in the H^{MY norm and equip this space with the classical H^{M) norm and inner 
product. 

The main linear portion of the equation is defined by^ 

-liT^T-vrdzzTy U={v,T,S)eD{A). 
-fis^S - vsdzzS) 

where we take: 

D{A) = {[/ = (v, T) e V(2) ■■ v^dzv + avv = 0, vtOzT + utT = 0, dzS = on T;, 

d^„T = d^„S = on TudzT = 9,5 = on Tfc}. 



^One sometimes also finds the more general definition [U, C/") := /v((v ■ v'd + ktTT'^ + K,$SS'')dAi with kt, Kg > fixed 
constants. These parameters Ky , Kg are useful for the coherence of physical dimensions and for (mathematical) coercivity. Since 
this is not needed here we take Ky = Kg = 1. Similar remarks also apply to the space V. 

^In [42] a slightly different non-selfadjoint operator also called A is given which corresponds to A + Ap (cf. (6.8)) in the 
present manuscript. 
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We observe that A satisfies the conditions given in Section 2.1. Note also that due to [49] (see also [42]) 
\AU\ = \U\h2. 

We next turn to the quadratically nonlinear terms appearing in (6.6). Noting that there is no momentum 
equation for w in (6.6) and in accordance with (6.6b) we define the diagnostic function: 

w{U) = w(v) = / V • vdz, U = (v, T, S) e V. 

J z 

Take, for [/,[/» (^ D{A): 

/ (v • V)v«\ / u;(v)a,v«\ 

Bi(U, [/») := Ph (v V)r» , B2{U, [/«) := Ph w;(v)5,T« (6.7) 

V(v-V)5V Vu;(v)5,W 

and let B{U, C/") := B{U, U^) + B{U, f/"). As in [42] one may show that B is well defined as an element in 
H for any !/,[/" G D{A). Furthermore B satisfies the conditions (2.3), (2.4), (2.5) relative to the definitions 
of A, D{A), V and H given here. For the second component of the pressure in (6.6a) we take 

/ -g/°(/3TVT + /35V5)dz 
ApU ^Ph\ ^ \, U eV. (6.8) 

and capture the Coriolis forcing in 

EU = Ph\ I , U C,H. (6.9) 

Finally we set 

Fu = Ph\ Ft I . (6.10) 





We may therefore define 

F[U) = ApU + EU + Fu (6.11) 

and observe that F : V ^ H and satisfies the requirement (2.8). Finally we define 

a{{v,T,S))^<j{U)=PH\ cjt{U) , UeH, (6.12) 

V '^s{U) I 

and assume either (2.10) or (2.11) for the consideration of Martingale or Pathwise solutions receptively. 

With the above definitions in place we may write (6.1) supplemented by the boundary conditions (6.2), 
(6.3), (6.4) in the abstract form (1.1) and conclude via Theorem 2.1: 

Theorem 6.1. Assume that F^, Ft and Fs are in Lf^^{[0,oo), L'^{A4)) and suppose that a{-) associated to 
o'v('); o't(')' '^s{') 'via (6.12) satisfies (2.10). Finally suppose that {vq,To, So) takes values in V and that 
/^o(') = lP((v,T', 5) e •) satisfies the moment condition (2.16) with q>8. Then: 

(i) There exists a local martingale solution of (6.1), the primitive equations of the ocean. 

(a) If we additionally assume (2.11) for a and allow of the relaxation of (2.16) to cover any q > 2 then 
there exists a unique maximal, pathwise solution of (6.1). 
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Remark 6.1. As noted in the introduction the Stochastic Navier-Stokes equations have been extensively 
studied. Initially these equations were considered with an additive noise. See [3] and later, [12, 35]. In this 
case, a classical transformation, allows one to treat u) € Q as a parameter in the problem. For nonlinear 
multiplicative noise, the Navier-Stokes equations were initially studied in the context of Martingale solutions. 
See, for example, [^8, 9, 8, 16, 37]. These works typically considered, from the PDE point of view, weak 
solutions evolving in time only in L^. . More recently, pathwise solutions (both local and global in time) for a 
multiplicative noise have been investigated. See [5, 2, 6, 36]. None of these works addressed the existence and 
uniqueness of a local pathwise solution in 3D evolving in H], . This case, whose deterministic analogue would 
be considered classical, was established only recently in [23]. As in the present work a key difficulty involves 
finding suitable compactness methods to pass to the limit. In [23] the necessary compactness is established 
by directly showing that the sequence of Galerkin solutions are Cauchy. For an application of this approach 
to the 2D Primitive equations see [24, 25]. In any case the work here provides an alternative proof of the 
results in [23]. The reader may readily check that the abstract framework developed in Section 2 applies to 
the Navier-Stokes equations on a bounded domain with Dirichlet boundary conditions and with H and V 
corresponding, approximately speaking, to L^ and H^ . 

7 The Passage to the Limit 

In this section we provide the details of the passage to the hmit, which is used in the proof of the existence 
of both martingale solutions and pathwise solutions. See Proposition 4.1 and Section 5.2 above. 

Proposition 7.1. Let Zk = {JJ'^'' ^W™''') be a sequence of X valued random elements mapping from a 
probability space (il, J^, P). We assume that 

(i) Zk converges almost surely to an element Z in the topology of X , i.e. 

{/"""^U inL^{[0,T],V)r\C{[0,T];V'), (7.1a) 

W'^'-^W inC{[0,T];iio). (7.1b) 

(ii) EachW"^^ is a cylindrical Wiener process relative to a filtration T]^'' that contains a{(W"^'' {s),U"^'' (s)); 
s<t). 

(Hi) Each pair {(}"■>' ,W"'') satisfies 

dU'"' + [AII'"' +9{\\U"'' -[/:''' II )B"'= ([/"'= )+i^"'=(^7"'')]di = cr"'=(^7"'=)dW^"^ 



if''"'{0)^Pn,Uo:^U(,'' 



(7.2) 



where we define U"'' by 

and assume, for some p > A, that 



^u:^ + Au:''^o ^7r(o) = ^7^'^ (7.3) 

at 



mUo\\P < oo. (7.4) 



Let S = (il, J^, {J-t}t>o^ Pi W), defining J-t as the completion of a{W{s), U{s); s < t). Then {S, U) is a global 
martingale solution of (2.22). Moreover if we define the stopping time: 

T:=inf |||C/-J7J| > 



1[\\U~U^\\>kY (7.5) 



where k is constant appearing in the definition of 9, (3.1), then (S,U,f) is a local martingale solution of 
(1.1). 
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The rest of this section is devoted to the proof which proceeds in stages. The first step is to estabhsh 
that the candidate sohition is in better spaces via (3.8) and weak compactness arguments. We next estabhsh 
the almost sure hmits of each terms arising in (7.2) against sufficiently 'smooth' test functions [/'* £ D{A). 
We then show, using (3.8) for [/"*= and the Vitali convergence theorem (see e.g. [20]) that each of the 
deterministic terms converges in L^(ilx [0, T]). The convergence of the stochastic terms in (7.2) are facilitated 
by Lemma 2.1. With these convergences in hand we make use of a variational argument (sc e.g. [23]) to 
finally conclude (2.22) for almost every time t and uj £ (l and with the equality understood in H. We pass to 
the limit for every t by establishing the improved continuity of U. See Subsection 7.3 below. This improved 
continuity justifies the definition of the stopping time r specified by (7.5). We therefore infer that for every 
t> 0, 

9{\\U -U4)B{U)ds^ B{U)ds. (7.6) 



In this manner we finally conclude that {S,U,t) is a local martingale solution of (1.1) and complete the 
proof. 

7.1 Improved Regularity of the Candidate Solution 

By applying the Banach - Alaoglu theorem with (3.8) for {/"*= in the case p = 2 wc infer the existence of 
elements tj e L^{n;L^{[0,T],D{A)) and U e L^{n:L°°{[0,T],V)) such that 

(/""^U mL'^{n;L^{[0,T],D{A)), (7.7) 

and 

(j^k_,*ij mL'^{n;L°°{[0,T],V). (7.8) 

On the other hand, due to (7.4), applied to Lemma 3.1, (3.8) we infer that for some q > 2, 

supE sup [£/"'= 1^, <c sup E sup ||C7"'=||« < oo. (7.9) 

k te[o,T] k te[o,T] 

Thus, with (7.1a), the Vitali convergence theorem implies that, 

U'^'^U \nL^{(l,L°°{0,T-V')). (7.10) 

Take Tl C [G,T] x Jl, measurable and [/« € D{A). By applying (7.7), (7.8), (7.10), we find 

E / XB.{U, U^)d3 = E f XBiU, U^)ds - E / xb^U, U^)ds, (7.11) 

Jo Jo Jo 

which means that U — U = U and wc conclude that 

U £ L^{n,L'^{[0,T],D{A)))nL'^{n,L'^{[0,T],V)). (7.12) 

Furthermore with (7.7) we have 

C/"" ^U in L'^in; L^{[0, T],D{A)). (7.13) 

7.2 Variational Equality for the Cutoff System 

Fix e D{A). Since, almost surely, [/"'= -^ U in L'^{[0,T], V) and noting that 

\ 2" 



supE 
fc 



\\ij"''fdt 

, 



<supcE sup ||t/"'=||'' < oo 
k \te[o,T] I 



24 



A. Debussche, N. Glatt-Holtz and R. Ternam 



we infer that [/"*= -> U in L'^{fl, i^([0, T], V)), by the Vitah convergence theorem. By thinning the sequence 
rik a necessary, we may also conclude that 



for almost every {t,uj) <E [0,T] x fl. 

The pointwise convergence in the linear term is direct: 



(7.14) 



{AiU''" -U),U^)ds 



<c\\U^\\ 



1/2 



[[/""-c/irds 



We conclude that for almost every {t, uj) E [0, T] x ft 



Jo 



(7.15) 



For B we estimate 

t 



(0(||c/"'= - ;/,"'= II )B"n^"'=) - ^(11^ - u4)B{u), u^)ds 







< 



< 



{eiWu"" - u^^lDiBiu^'") ^ B{u)),Pn,u*)di 



{{eiWu''" ^ ui'"]]) - e{\\u - u4))Biu),u^)ds 

{eiwu^^" -U:-\\)B{U),Q,,,U^)ds 



mwu"" - u:^ ID - em- u4))b{u), u^)ds 

We address the elements on the right hand side in reverse order. Due to (7.14) we have, for almost every 

(t,w)e [o,r] xn-.^ 

We therefore infer, 

d{\\u"-u:\\)^e{\\u-Lf4), 

almost everywhere on [0,T] x ft. By assumptions (2.3), (2.4), 



((0(11^7"'= - f/r^ll) - 0{\\U- U4))B{U), [/«) dt < c\AU^\\\Uf, 



and since. 



E 



T nt 



"'0 



m\\U''^ - f/r II) - 9{\\U - U4))B{U), [/«) 



dsdt 



< cE 







U">'~U':'\\)-0{\\U-U4))B{U),U^) 



dt, 



the Lebesgue dominated convergence theorem therefore implies: 

E / J'^'-dt^Q. 



^Since C/J -> Uo in L'^{n; V), it follows that E (sup^gjo ^j \\IJ^ - U,\\A -> 0. 
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Thinning the sequence, if necessary we conclude that 

JJ^" ^ a.e. (t,uj) e [0,T] x Cl. (7.16) 

We estimate the second term, Jj'*" according to 

.72" < \Qn,AU*\ WUfds^O, (7.17) 



JO 

for almost every {t,Lu). Finally for J"*", the bilinearity of B implies that 

Again by the assumptions (2.3), (2.4) we infer, 

Ji"^ <c\AU^\ f (||i7"^|| + WUIDWU''" - U\\ds 
Jo 

<c|AC/«| [ /" (||C/"^f + ||t/|p)rfsj If WW'^-UW'ds] . 

Thus, with assumption (7.1a), we infer that 

J"" -^ for almost ah (i,a;). (7.18) 

Combining (7.18), (7.17), (7.16) we conclude that, for almost every (t,cj), 

(e'di?/"'" -t/r'=ll)B"(c/"),[7'')ds^ / {e{\\u-u4)B{u),u^)ds. (7.19) 

"'0 

For the remaining deterministic terms we estimate 



<c|(7«| / \FiU^^)-FiU)\ds + \Qn^U^ \F{U)\ds:=J2'+J^''. (7.20) 



(F"'=(C/"'=)-F(t7),t/")ds 

Due to (7.14) and the continuity assumed for F, (2.7), we infer that for almost every (w, t), 

iFid"") - F{U)\ -^ 0. 

On the other hand by (2.7), 

iFiU""-) - F{U)\ < c(l + ||[7"^|| + ||C/||), 

and we infer that, 

supE /" \F{U'"') - F{U)\^dt < csupE f (1 + Ht/^'-f + ||t/|p)dt < oo. 

k Jo k Jo 

In consequence {|F([/"*^) — F{U)\}k>o is uniformly intcgrable over Ct x [0,r]. By applying the Vitali Con- 
vergence Theorem we have E/„ |i^"fc({7"'!) _ F{U)\dt — > 0. Thinning the sequence further if needed, we 
infer that almost everywhere in f2. 

f \F''-{U'"-)~F{U)\dt< [ \F'"'{U'"')-F{U)\dt^O. 
Jo Jo 

in order to finally conclude that 

Jl'^^O a.e.{uj,t). (7.21) 
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Turning to the second term J^'' , we see, again as a consequence of the assumption (2.7), that /„ \F{U)\ds < 
c/p (1 + ||f7||)(is < oo and so 



J^-^O a.e.{io,t). 



In conclusion, by (7.21), (7.22) we finaUy have 



f {F''-{U''-),U^)ds^ [ {F{U),U^)ds 
Jo Jo 



(7.22) 



(7.23) 



for almost every (w,t) € fi x [0,T]. 

We next establish the convergences to the deterministic terms in (2.22) in the space L'^{il x [0,r]), 
1 < g < 2. Notice that due to (2.4), (2.7), 



E 



(AC/"" +9{\\U"'' - U:^''\\)B''{U'"') + F"" {&""), U^)ds 



dt 



<cE / liAU''" + eiWU''" - Ul"'\\)B'\U'"') + F'\U'"'),U^)\^ds 



<c|AC/«pE / (||;7"'=f + wu'^-w^ + l)ds. 
Jo 

Thus, for every q £ [1, 2), 

f (At/"" +6l(||t/"'= -C7,""||)B"'=({/"'=) + F"'=([7"'=),[/»)ds 

'0 J k>0 

is uniformly integrable in L'^{i} x [0,T]). 
Combining this with (7.15), (7.19), (7.23) we conclude that for every q S [1,2): 

(A?/"" +6i(|l;7"'= - C/,"'"||)B"'=(;7"'') + F"''-(C/""),[/«)ds 

{AU + ei\\U - U,\\)B{U) + F{U), U^)ds, 



in L''{[0,T] X n). 

The stochastic terms are handled differently. Using (2.2) and (2.10) wc estimate 

||a"" (?/"") -a(f7)|U,(u,i,)<||a(t/"")-a(f/)|U,(u,H) + ||Q„.a(t/)||i,(H,ff) 

<i|a(L/"'=) - <7(;7)|U,(u,H) + -^|k(f7)i|L.(u,v) 

<||a(C7"") - a([/)|U,(u,H) + -^(1 + l|t/||). 

Thus, due to (7.14) and the assumed continuity of a (see (2.10)) we conclude that 

||a"'=(f7"'=)-a(f7)|U,(u,H)^0, 
for almost every (w, t) G il x [0, T]. On the other hand, we observe that 



supl 



"^ \K'{U''')\\Un.H)dA < csupE ( ^(1 + \\U-^f)ds] , 

J nk \Jo I 



(7.24) 
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where we have again made use of the subhnear condition (2.10). We therefore infer that |icr"'"(C/"'°)||L,(u,ff) 
is uniformly integrable in LP{il, x [0, T]) for any p £ [I, 4). With the Vitah convergence theorem we infer, for 
aU suchp e [1,4), 

<j"''iU"'')^aiU) inLP(fi;LP([0,T],i2(it,i7))). (7.25) 

In particular (7.25) implies the convergence in probability of cr"'= ([/"*=) in i^([0, T], L2{il, H))). Thus, along 
with the assumption (7.1b), we apply Lemma 2.15 and infer that 

cr"" ({/"'-- )diy„, -> / a{U)dW, (7.26) 

"'0 

in probability L^([0, T], H)). Another application of the Vitali convergence theorem using estimates involving 
(2.12), (7.25) shows that the convergence in (7.26) occurs moreover in L^(J7; L^([0,T], iJ)). 

With the above details in hand we now establish (2.22) in a variational sense. Fix any C/" e D{A), 
7?. C O X [0,r] measurable. Using (7.13) and then (7.24) and (7.26) wc observe that 

E / X'R{U,U^)dt^ lim E / X7^(^7"^ [/»)di 
Jo '-■^°° Jo 

= lim E / xn{Ui;\U^)dt 
^^°° Jo 

- lim E / xn\ I (A^7"^t/«)ds ) dt 
fe->oo Jo V^o / 

- lim E / xn\ I {OiWU""- - f/:'^ ||)B"^ ([/""), t/»)ds) dt 

fe^oo Jo \J0 / 

- lim E / xn\ I {F"'{U"''),U^)ds] dt 
>''^°° Jo \Jo J 

+ lim E / xni I {(T'"'{U'"'),U^)dW'"'] dt 
f"^^ Jo \Jo J 

( Xn({Uo,U^)- I {AU + 9{\\U-U4)B{U) + F{U),U^)ds]dt 

+ eJ xn(j {a{U),U^)dw\dt. 

Since this equality holds over all such TZ we may conclude that for almost every (w, i) G J7 x [0, T] and every 
G D{A) that, 

([/(t),L/«)+ / {AU + 0{\\U -U4)B{U)+F{U),U^)ds = {UQ,U^)+ f {a{U),U^)dW. (7.27) 

Jo Jo 

Moreover, due to (7.12) established above, it follows by density that (7.27) holds also over [/" e iJ and hence 
(2.22) in the analogous sense to (2.19). 

7.3 Improved Regularity In Time 

With (7.27) and (7.12) in hand it remains only to establish better continuity, in time, for U. More precisely, 
we must show that U G C'([0, T]; V) a.s. Of course, such a condition is needed in order to justify the definition 
(7.5). 

To this end wc define 

dZ + AZ = a{U)dW, Z{0) = Uo. (7.28) 
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Observe that since a{U) £ L^{n,L^{[0,T],L2{ii,V))) we have 

Z e L\n, C([0, T],V)) n L^{n, L\[0, T];D{A))). (7.29) 

Now take U ^U - Z. Subtracting (7.28) from (2.22) we find 

j^U + AU + e{\\U + Z-U4)B{U + Z)+FiU + Z) = 0, ^^^^^ 

t/(0) ^ Uq. 

Due to (7.29), (7.12), we infer that U e L'^{n, L^{0,T; D{A)) n L°°{[0,T],V) and hence that, 

AU, 9{\\U + Z-U4)B{U + Z), F{U + Z) e L^in,L^{[0,T],H). (7.31) 

We conclude with (7.30) that 

^A^^^U eL^in;L^iO,T;V')), A^/^U e L\n;L^{0,T;V)). (7.32) 

By applying [47, Chapter 3, Lemma 1.2] we infer that A^/'^U G C{[0,T],H) so that, with (7.29), we deduce 
that 

il e C{[0,T];V), a.s. (7.33) 

With (7.27), (7.12), and (7.33) we finally conclude that {S,U) is a global Martingale solution of (2.22). 
Furthermore, having justified (7.5) and applying (7.6) to (7.27) we have that {S, U,t) is a local Martingale 
solution of (1.1). The proof of Proposition 7.1 is therefore complete. 

8 Appendix: Proof of the Convergence Theorem 

In this final section we provide a proof of Lemma 2.1. Convergence results similar to Lemma 2.1 have 
appeared in previous works (see e.g. [1], [28]). However, to the best of our knowledge, no one up to the 
present has provided a detailed proof. Note that in the present work Lemma 2.1 is an important technical 
tool for the passage to the limit, as detailed above in Section 7. 

To simplify the exposition, we begin by introducing the notations: 

For the truncations we set 



T^ '■— 2^ ^h 1 Jn '■— 1-^ - 1-N ^ 1-N '■— 2^ Yk, Ji^:—I^-I\ 



n Tn 

N- 



N>k>0 N>k>0 

With these notations we now split 

\^" -^\l-^{[o,t],x) < |2^" -^n\l^{[o,t],x) + \^N ~^n\l'^{[Q.T],X) + l^N -^\l^([o,t],x) 
and observe that the proof of Lemma 2.1 is complete once we establish that 

For every e > 0, lim sup P {\Jn\lV[o,t].x) > e) = 0, 

lim |y," — Yk\L2riQ rpi j^s — in Probability, for each fixed k, (8-1) 

lim \JN\L^no,T],x) = in Probability. 
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To establish each of the convergences in (8.1) we make extensive use of the following martingale inequality 
(see e.g. [21]) 



FdW 



X 



dt>c\ < — 



\F\ 



L2iii,X) 



dt > K 



■2) 



Here c, k may be any positive constants and F any Tt predictable element in L'^{[0,T];L2{ii,X)). For the 
first item in (8.1) we apply (8.2) and observe that for any e,^ > 



5 
< - 
- 3 



With this estimate, the assumptions on G and (2.14b) we infer the first item in (8.1). The third item in 
(8.1) is established in similar manner via an application of (8.2). 

It remains to address the second item in (8.1). In order to treat these terms we introduce the functional: 



T^piP) = - [ cxp (-- — -) F{s)ds F e L\[0,T],X),p> 0. 
P Jo \ PI 



(8.3) 



Using this functional and then integrating by parts we estimate 



\yH-yAx = 



< 



< 



GldWl} - / GkdWk 

JO 

t 



X 



{Gl n,{Gl))dW^ 



f(Gl^n,{Gl))dWl^ 
Jo 



X 



X 



{TZpiGk) - Gk)dWk 



(UpiGk) - Gk)dWk 



X 



X 



UpiGDdWj:- / 'Rp{Gk)dWk 



X 



(8.4) 



\'Rp{Gk)W^~np{Gl)W]:\x 



{np{Gu)Wk - 'Rp{Gl)W]:)ds 



1 

+ - 
X P 



{GkWk ~ GlW]})ds 



We now proceed to treat each of the term on the right hand side of (8.4). Fix e, J > 0. For the first term 
in (8.4) wc apply (8.2) and estimate 



{Gl-np{Gl))dW^ 



L2([0,T];X) 



>€\<6 + : 



5j^ 
T 



<S- 



\G]^ - Gk\Ut > ^— I +] 



ST 



<S 



ij^\G^^^np{Gl)\\dt> 

(£\Gk-npiG,)\Ut>^^] 

(£\np{G,)^TZp{Girxdt>^^] 



(8.5) 



\G",~Gk\Ut> 



With (8.2) wc also find that 



(TZpiGk) - Gk)dWk 



L2([0,T];X) 



>e] <6 + : 



\Gk-npiGk)\'Ut> 



5j^ 
T 



(8.6) 
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The last three items are treated differently 



' T ,2\ ( i-T 



<P n \ll,{Gu)Wl -n,(Gl)W^\\dt > M +P n \n,{G^)W^-Tl,{G^)Wi:\\di > '-] (8.7) 

<P ( sup \W^\' f \Gu - GlWdt > ^ ) + P ( sup |M^, - W]:\^ f \Gu\\dt > f ) . 
\te[o,T] Jo 4 / \te[o,T] Jo ^ J 



Similar estimates lead to 



- f {TZ,{Gk)Wk^TZp{Gl)W;!)ds 
P Jo 



> e 



L^i[0,T];X) 



sup \Wl\^ \ \Gu GlWdt > lA + p sup |vi/, - Wl\^ \ \Gu\\dt > 1^ 
tG[0,T] Jo 4J J \telo,T] Jo 4J 

The final term in (8.4) yields to an identical estimate. 

Collecting the estimates (8.5), (8.6), (8.7), (8.8) wc infer that 

p(|r,"-yfeU.([o,T],x)>5e) 



<2S + 2P ( / \Gl - Gfcli-rff > _ ) + 2P ( / \Gk- np{Gk)\ldt > 



(8.9) 
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sup \W^\' f \Gk - GlWdt > ^1 +3P ( sup \Wk - Wn' f \Gk\ldt > ^ 

^tG[0,T] Jo 4J J \ti£[0,T] Jo 4i 



Since S,e > are arbitrary and given basic properties of the functional (8.3) along with (2.14) we may now 
infer the second item of (8.1) from (8.9). This completes the proof of Lemma 2.1. 
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